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Abstract 

By using the scaling method we derive the virial theorem for the relativistic mean field 
model of nuclei treated in the Thomas-Fermi approach. The Thomas-Fermi solutions 
statisfy the stability condition against scaling. We apply the formalism to study the 
excitation energy of the breathing mode in finite nuclei with several relativistic parameter 
sets of common use. 
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The relativistic mean field (RMF) treatment of Quantum Hadrodynamics [[I], [2j has proven 
to be very useful for describing different properties of nuclei along the periodic table. The 
simplest model, the linear o — uj model of Walecka 0, explains the nuclear force in terms of the 
exchange of o and uj mesons. It is known that the value of the nuclear matter incompressibility 
is unreasonably high in this linear model (K nm ~ 550 MeV). The problem can be cured by 
introducing cubic and quartic self-interactions of the a meson Q , and the model can be refined 
by adding an isovector p meson. Current non- linear parameter sets, such as the NL3 set |J, 
give ground-state energies and densities in excellent agreement with the experimental data, 
not only for magic nuclei but also for deformed nuclei and for nuclei far from the stability line. 

The scaling method has been often employed to derive the virial theorem in the non- 
relativistic framework ||, e.g., for nuclear effective interactions such as the Skyrme force [0]. 
It has also been applied in calculations of nuclear collective excitations like the breathing mode 
(isoscalar giant monopole resonance) J/J. Relativistic generalizations of the virial theorem 
obtained by use of the scaling method exist for particles in external potentials || [JJ. In the 
RMF model of nuclei the mean field potentials are generated self-consistently. Owing to the 
meson-exchange nature of the relativistic model one has to deal with finite range forces, which 
renders the scaling more involved than for zero-range Skyrme forces. Moreover, in contrast to 
the non-relativistic situation, there exist two different densities, namely the baryon and the 
scalar density, in accordance with the fact that one has two types of fields, the vector and the 
scalar field. 

In this letter we shall make use of the principle of scale invariance to obtain the virial 
theorem for the RMF theory by working in the Thomas-Fermi approximation. We shall 
include non-linear self-couplings of the scalar field and shall deal with spherical finite nuclei. 
The second derivative of the scaled energy with respect to the scaling parameter, the so-called 
restoring force, turns out to be positive (stability condition) in the Thomas-Fermi calculations. 
Thus we are able to apply the method to compute the excitation energy of the isoscalar giant 
monopole resonance in finite nuclei with realistic parameter sets of the relativistic model. 
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The meson field equations of the non-linear a — uj model are 
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Here p = p p + p n is the baryon density, p 3 = \{p p — p n ) is the isovector density, and p s is 
the scalar density. The meson fields <f), V and R are associated with the a, uj and p mesons, 
respectively, and A is the Coulomb field. It is understood that the densities and fields are local 
quantities that depend on position, even if we do not make it explicit. Units are h = c = 1. 

Taking into account the above field equations, the relativistic energy density of a finite 
nucleus in Thomas-Fermi approximation can be written as [|1], [4| 



H = £+ l -g s 
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in terms of the nucleon energy density 
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is the scalar density and m* = m — g s <p is the nucleon effective mass. For each kind of 
nucleon (q = n, p) the local Fermi momentum k Fq is defined by kp q = (37r 2 p g ) 1//3 , while 



ki n + m . 



e Fq - \/ ^Fq 

The virial theorem relates the kinetic and potential energy components of the energy in 
certain circumstances. This theorem results from homogeneity properties of the kinetic and 
potential energy components of 7i with respect to a scaling transformation that preserves the 
normalization. One such normalized scaled version of the baryon density is 



p x (r) = A 3 p(Ar), 



(8) 



where A is an arbitrary scaling parameter. Accordingly, the local Fermi momentum changes 
as 

k Fq x(r) = [3tt 2 P(? a(7-)] 1/3 = Xk Fq (Xr). (9) 

The meson fields and the Coulomb field are also modified by the scaling due to the self- 
consistent equations (|I])-(f4|), which will relate the scaled fields to the scaled densities. Un- 
fortunately, the meson fields do not scale according to simple power laws of A because of the 
finite-range character of the meson interactions. This is most apparent for the scalar field (f>, 
since the scalar density in the source term of Eq. ([!]) transforms not only due to the scaling 
of kp q but also of itself (or m*), see Eq. (0) for p s . For reasons that will become clear 
immediately, we shall write the scaled effective mass m* x (r) = m — g s <p\(r) in the form 



m x (r) = Am*(Ar). 



(10) 



The quantity fh* carries an implicit dependence on A apart from the parametric dependence 
on Ar. 

On account of Eqs. @ and (|10|) the scaled form of £ reads £\{r) = X A £[kF q (Xr), m*(Xr)] = 
A 4 £(Ar), while the scaled scalar density reads p s \{r) = X 3 p s [k Fq (Xr) , m*(Ar)] = A 3 p s (Ar). The 
tilded quantities £ and p s are given by Eqs. (H) and (0) after replacing m* by in* . Note the 



usefulness of fllCf ) to be able to put the transformed densities £\ and p s \ into the above compact 
form. This way, for the scaled total energy density 7i\ we obtain 



X" 
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with the definition gsPt s = g s p s — b4>\/X 3 — c0|/A 3 . 

The scaled energy is stationary for A = 1 (which leads to the virial theorem): 
d rd{Xr) 
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Here we have used dS/dX = p s dm*/d\ (as p s = dS/dm*) and, from the definition of m*, 
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Let us exemplify the calculation of the derivatives of the scaled fields with respect to A 
with the omega field V\. It fulfils the scaled Klein-Gordon equation (A u — m 2 /X 2 )V\{u) = 
—Xg v p(u), where we have used Eq. (|§D for p x and have switched to the coordinate u = Xr. 
On differenciating this equation with respect to A we have 
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If one now sets A = 1 the solution of this equation provides dV\/ dX\\ = \. Nevertheless, for our 
purposes it is more useful to multiply both sides of ( |i~4"D by V\, integrate over the space and 
then use Green's identity on the left hand side. This way it is straightforward to get 
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which at A = 1 is just one of the contributions we need in Eq. (p!z|) . Analogous results 
are found for the rho and Coulomb fields (with a zero mass for the latter). In the case of 
the scalar field additional terms appear due to the fact that the scalar density itself is a 
function of the scalar field. Following the same steps as above, from the scaled field equation 
(A u — ml/X 2 )<f)\{u) = —Xg s p~l s (u) one easily arrives at 
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From substitution of Eqs. flT5|) and (O) (and of the corresponding results for the rho and 
Coulomb fields) into Eq. (O) the virial theorem for the non-linear a — uj model becomes 
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Actually, introducing the kinetic energy density r we have £ — m*p s = r + mp — mp s , which 
makes more obvious the kinetic energy component in the virial theorem. Using Eq. to 
eliminate S from the expression of the relativistic energy density 7i, the RMF energy of a 
nucleus takes the remarkably simple form 

1 



J dr[H — mp] = J dr 
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where we have subtracted the nucleon rest mass contribution. This expression shows very 
clearly the relativistic mechanism for nuclear binding. It stems from the cancellation between 
the scalar and vector potentials and from the difference between the scalar and the baryon 
density (i.e., from the small components of the wave functions). Equations fljjp and (|19D are 
satisfied not only by the Thomas-Fermi solutions, but also by the ground-state densities and 
meson fields obtained from a quantal Hartree calculation. Of course, the energy stationarity 
condition against dilation of the RMF problem must be fulfilled by any approximation scheme 
utilized to solve it. 

As a further application of the method we shall use it in calculations of the isoscalar giant 
monopole resonance (ISGMR). It is customary to write the excitation energy of the ISGMR 
as 

[FT. 

(20) 
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where Cm and .Bm are called, respectively, the restoring force (or incompressibility of the finite 
nucleus) and the mass parameter of the monopole vibration. To study Em in the RMF the 



authors of Refs. flT0| , |TT| resorted to a local Lorentz boost and the scaling method. Following 
these works one has 
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where the scaling parameter A now plays the role of the collective coordinate of the monopole 
vibration, and 

B M = jjdrr 2 H(r), (22) 
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with A being the mass number of the nucleus. The investigations of Refs. |T0| , [TT| were 



restricted to the linear a — to model, either for nuclear matter or for symmetric and 



uncharged finite nuclei (with the densities solved in Thomas-Fermi approximation) fll |. It is 
well known that the surface properties of nuclei cannot be described within the linear model, 
and hence nor can the overall properties of nuclei. 

To compute Cm it is easiest to replace the relations (|1"5| ) and fllTD into the expression 
of d[J drH\(r)]/d\ and derive again with respect to A. After some algebra we obtain the 
restoring force as 



C 



M 



.4 



dr 



— m- 



dp s 
dX 



+ 3 (m s 2 2 + hcf - m\V 2 - m 2 p R 2 ^ 



(2m 2 + 60 2 )^ > ^ 2 

OA 



+ 2miV^ + 2mlR 9Rx 



dX 



J A=l 



where 



dp s 
dX 



dp s dm* 



A=l 



dm* dX 



dp s 



A=l 



dm* 



m* + g s 



dX 



(23) 



(24) 



J A=l 



The derivatives of the scaled meson fields with respect to A are computed by solving Eqs. (|i~4T) 
and flT6| ) at A = 1. We have found Cm to be positive for all of the (linear and non-linear) 
parameter sets we have tested in the Thomas-Fermi calculations. A large part of the final 
value of Cm (usually far more than a half) is due to the contribution of the term dp s /dX\\=i- 
In Table 1 we display the calculated Thomas-Fermi excitation energies of the ISGMR, 
together with the empirical estimate E M ~ 80/ A 1 / 3 MeV [0, for 40 Ca, 90 Zr, 116 Sn, 144 Sm 
and 208 Pb. Recent experimental data on the centroid energy of the ISGMR are available for 
these nuclei [jHJ. We have employed the non-linear parameter sets NL-Z2 (K nm = 172 MeV) 
rj, NLl (K nm = 212 MeV) |T5|, NL3 (K nm = 272 MeV) §, NL-SH (K nm = 355 MeV) 
16| and NL2 (K nm = 399 MeV) ||17| . These parameter sets have been determined by least- 



squares fits to ground-state properties of a few spherical nuclei and are of common use in 
RMF calculations. From the table one can see that the smaller the mass number, the larger is 
the monopole energy. The energy of the ISGMR increases with increasing K nm in the various 
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parameter sets. For example, the monopole energy in 208 Pb is 12.3 MeV for NL-Z2 while it is 
18.1 MeV for NL2. The dependence on K nm is roughly linear for each nucleus. 



In assuming nuclear matter within a certain volume the authors of Ref. QTUJ] estimated the 
monopole excitation energy of a finite nucleus as 



Em — \ T~2U ~T — t7 — V' ( 2 ^) 

V {r 2 ){e F nm + g v V nm ) 

with (r 2 ) = |i? 2 and R = 1.2v4 1//3 fm. They evaluated ( |25|) for the linear model of Walecka 
and found Em = 160/A 1 / 3 MeV, which has the correct dependence on the mass number but 
is twice as large as the empirical value. From Eq. ( pq) one finds Em = 92, 102, 115, 132 
and 140/A 1 / 3 MeV for the non-linear sets NL-Z2, NL1, NL3, NL-SH and NL2, respectively. 
Comparing with Table 1, the finite size effects reduce the prediction obtained from nuclear 
matter by a factor ranging from ~ 1.4 in 40 Ca to ~ 1.3 in 208 Pb, rather independently of the 
parameter set. 

The ISGMR has been studied in the time-dependent RMF (TDRMF) theory by Vretenar 
et al. ||18|| . We include in Table 1 their results for the energy of the main peaks that appear in 
the monopole strength distributions of 90 Zr and 208 Pb. Our scaling results compare very well 
in the case of 208 Pb for all parameter sets, but give somewhat larger excitation energies for 
90 Zr. It should be mentioned that the Fourier spectrum of 90 Zr in the TDRMF calculation is 
considerably fragmented (specially for the sets with higher K nm ) and then the determination 
of the centroid energy is more uncertain [|TS| . Very recently, it has been demonstrated that the 
relativistic random phase approximation (RRPA) is equivalent to the small amplitude limit 
of the TDRMF theory in the no-sea approximation, when pairs formed from the empty Dirac 



sea states and the occupied Fermi sea states are included in the RRPA [T9| . 

Microscopic calculations of ISGMR energies in nuclei are a valuable source of information 
on the nuclear compression modulus K nm [^, [HJ, which is an important ingredient not only 
for finite nuclei but also for heavy ion collisions, supernovae and neutron stars. A further 
inspection of Table 1 shows that the empirical law Em ~ 80 /A 1 ^ 3 MeV lies between the 
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predictions of the NL1 and NL3 sets, as expected from the reasonable value of K nm in these 
interactions. On the contrary, K nm is too high in NL-SH and NL2 and we see that these 
sets overestimate the empirical curve and the experimental data for all nuclei of Table 1. 
No RMF parameter set seems capable of reproducing the mass-number dependence of the 
experimental data over the whole analyzed region, particularly in the lighter nuclei. One 
should note, however, that our calculation provides a prediction for the mean value or centroid 
of the excitation energy of the resonance. To establish a link between K nm and the measured 
energies the most favourable situation is then met in heavy nuclei, where the experimental 
strength is less fragmented than in medium and light nuclei If we only take into account 
the data of 144 Sm and 208 Pb, our results of Table 1 suggest that K nm of a RMF interaction 



should belong to the range 225-255 MeV. (If we disregard the set NL-Z2, as in Refs. |T_8j, |TJ 



the range is 230-260 MeV.) From their TDRMF and RRPA calculations the authors of Refs. 



TB| , [TJ| conclude that the value of K nm should be close to 250-270 MeV. Non-relativistic 



Hartree-Fock plus RPA analyses using Skyrme and Gogny interactions determine K nm to be 



215 ± 15 MeV gg, |JJ, thus lower than in the RMF model. 

We have derived the virial theorem for the relativistic nuclear mean field model on the 
basis of the scaling method and the Thomas-Fermi approximation. In this approach we 
have calculated for realistic parameter sets of the RMF theory the breathing-mode energy 
of finite nuclei fully self-consistently (i.e., we did not use a leptodermous expansion of the 
finite nucleus incompressibility as in some previous studies with the scaling method |25|). The 



present calculations extend earlier work performed with the linear a — u model [10, 11 



The excitation energies of the monopole oscillation turn out to be in good agreement 
with the outcome of dynamical time-dependent RMF calculations. It has been shown very 
recently that the relativisitc RPA, with the inclusion of Dirac sea states, amounts to the limit 
of small amplitude oscillations of the TDRMF theory flT9"| . From the present Thomas-Fermi 
analysis one can thus conclude that, similarly to the non-relativistic case, also in the relativistic 
framework the excitation energies obtained with the scaling method simulate the results of 
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the random phase approximation. 
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Table 1: Excitation energy of the monopole state (in MeV) obtained in the scaling approach 
by using various relativistic parameter sets (in order of increasing value of the compression 
modulus K-n m ). The energies of the main peaks found in the time-dependent RMF calculations 



of Ref. |L8| are also shown for Zr and Pb. The experimental centroid energies are from 



Ref. 13 







NL-Z2 


NL1 


NL3 


NL-SH 


NL2 


80A- 1 ' 3 


Exp. 


40 Ca 20.5 


21.2 


23.5 


26.6 


29.5 


23.4 


19.2 ±0.4 


90 Zr 16.4 


17.2 


19.2 


21.9 


24.0 


17.9 


17.9 ±0.2 




18 




15.7 


~ 18 










116 Sn 15.1 


15.9 


17.7 


20.3 


22.3 


16.4 


16.1 ±0.1 


144 


Sm 14.1 


14.9 


16.6 


19.0 


20.8 


15.3 


15.4 ±0.3 


208 


Pb 12.3 


12.9 


14.5 


16.6 


18.1 


13.5 


14.2 ±0.3 
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17.8 
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